Charge pumping in monolayer graphene driven by a series of time-periodic potentials 
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We applied the Floquet scattering-matrix formalism to studying the electronic transport proper- 
ties in a mesoscopic Dirac system. Using the method, we investigate theoretically quantum pumping 
driven by a series of time-periodic potentials in graphene monolayer both in the adiabatic and non- 
adiabatic regimes. Our numerical results demonstrate that adding harmonic modulated potentials 
can break the time reversal symmetry when no voltage bias is applied to the graphene monolayer. 
Thus, when the system is pumped with proper dynamic parameters, these scatterers can produce a 
nonzero dc pumped current. We also find that the transmission is anisotropic as the incident angle 
is changed. 
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I. INTRODUCTION 

Quantum transport in periodically driven mesoscopic 
systems is currently of great interest because of an in- 
creasing number of applications. [ll-fl2| Many research 
studies have been carried out on the mechanism of 
scattering by time-dependent potentials in recent years. 
There exists a good theoretical understanding of the tun- 
nelling events in the presence of an alternating (AC) field, 
i.e., the Floquet theory, [l3l - [T6| which takes the photon- 
assisted transport into account and allows us to con- 
vert the solution of a time-periodic Schrodinger equation 
into a time-independent eigenvalue problem. For sys- 
tems with time-periodic potentials, inelastic tunnelling 
occurs when electrons at appropriate incident Fermi en- 
ergies undergo transitions between the central band to 
several sidebands by means of photon emission or ab- 
sorption, referred to as photon-assisted tunneling (PAT) . 
The central band just corresponds to particles with the 
incident Fermi energy while the sidebands correspond to 
particles which have gained or lost one or more mod- 
ulation energy quantum. Especially, if the potential is 
harmonic in time, the PAT results in exchange of energy 
with electrons in units of modulation quantum fiw, with lj 
being the modulation frequency. An important applica- 
tion of the PAT is quantum charge pump, which received 
considerable theoretical investigation and had been ob- 
served experimentally. 0-Efl E GJ] However most 
previous studies in this realm have been focused on the 
Schrodinger system. 

In Ref . fl(| . the Floquet scattering theory was devel- 
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oped for quantum pumping in mesoscopic semiconduc- 
tors by Moskalets and Biittiker. This approach permits 
description of both adiabatic and non-adiabatic regimes 
of pumping. It is the purpose of this work to applied the 
Floquet scattering-matrix method to calculating quan- 
tum charge pumping in the Dirac system. To this end, we 
extend the approach of Ref. to accommodate a novel 
two dimensional hexagonal carbon material, graphene. 
Recently, the study of quantum pumps in Dirac system 
has been carried out in one dimensional Luttinger li quid 
model [ItJ or weak pumps in the adiabatic regime [l8l.ll9l] 
and nonadiabatic regime (2p| . [2lj in graphene by using 
the Floquet scattering theory. Graphene, a single layer 
of carbon atoms arranged in a hexagonal lattice, exhibits 
abundant new physics and potential applications. (22r - 
[26| Quantum transport properties in graphene have at- 
tracted increasing attentions both for the fundamental 
physics and potential applications in carbon-based nano- 
electronics devices. The novel properties arise from the 
unique linear dispersion and the chiral nature of its car- 
riers near the Dirac points. At the K {K ) point of the 
Brillouin zone, the energy spectrum exhibits a linear dis- 
persion that can be well described by the massless Dirac 
equation with an effective speed of vf « 10 6 ms _1 . (2(| 
These quasiparticles, called massless Dirac fermions, are 
quite different from the standard electrons that we en- 
counter in conventional two dimensional electron gases 
in semiconductor heterostructures. Graphene monolayer 
is a suitable testbed to examine the photon-assisted elec- 
tron tunneling |27l l28j and quantum pumps in the Dirac 
system. In addition, the high electron mobility and long 
phase coherence present in graphene could in principle 
enhance these interference effects (photon-assisted tun- 
neling and quantum pumps) , thereby making possible its 
observation with the present technology. We carried out 
a study of quantum pumps in both adiabatic and non- 
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adiabatic regimes in graphene monolayer reported here, 
which helps to obtain a clear physical picture about time- 
dependent tunneling and will be interesting for the po- 
tential application of carbon-based electronic devices. 

In this work, we study theoretically the charge pump 
by a series of time-periodic potentials in monolayer 
graphene with no external bias. We also find exchange 
of photons between the oscillating barrier and electrons. 
Such processes give rise to transitions from the central 
band to sidebands, all of which can convey the elec- 
trons and contribute to the pumped current. The charge 
pumping is a consequence of the interference of differ- 
ent propagating modes of which the differences in phase 
break the time reversal symmetry. 

This paper is organized as follows. In Sec. II, 
we present the theoretical model and the calculation 
method. We show the numerical results and discussions 
in Sec. III. Finally, we give the conclusions in Sec. IV. 



II. THEORY 



A. Floquet scattering model in graphene 



ihftfd) - VI cos(wt + ^)f(t) = (E- E F )f{t), (5) 
where E is a constant. Integrating Eq. (JSJ) gives: 

f(t) = e~' l(E ~ E ^ t ' h eM-T /*t?C0B(wt +(p i )dt) 
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where we have taken /(0) = 1, and J„ is the Bessel func- 
tion of the first kind in order n. Note that f(t) = f(t+T), 
Eq. (IB]) requires that E = Ep + mhij, where m is an inte- 
gral index of the sidebands. The translational invariancc 
along the y direction gives rise to the conservation of k y , 
and thus the solutions can be written as g(r) = g(x)e lkyV . 
The equation for g(r) has a solution: 
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We consider electrons transmitting through a series of 
modulated pumping potentials in a monolayer graphene 
sheet. The charge pumping potentials can be realized by 
top metallic gates [29| or by using the electric field of a 
surface acoustic wave. [3(| We approximate the scatter- 
ing potentials by a sequence of time-periodic rectangular 
barriers. The barriers are infinite along the y-direction 
and homogeneous in each region as illustrated in Fig. [1] 
The potentials are given by: 

V l (r,t) = + VI ■ cos(wi + (4 < x < x^ 1 ), (1) 

where the superscript i indicates the i-th region (see 
Fig. [T]). The height of the barrier is oscillating sinu- 
soidally around V l with amplitude VI, frequency lo and 
phase ip l . Quasiparticles in graphene are formally de- 
scribed by the Dirac-like Hamiltonian: 



H = hv F (T -k + V\r,t). 



(2) 



where k is the momentum, a{i = x, y) is the pseudospin 
Pauli matrix, and the k- independent Fermi velocity vf 
plays the role of the speed of light. 

The Floquet theorem asserts that the solution of the 
eigen-energy problem in our system has the form: 



(3) 



where Ep is the Floquet eigen-energy and <f)(r,t) = 
4>{r, t + T) is a periodic function, with period T = 27r/w. 
Since the time and position dependence of V l (r,t) are 
separated in our model, <fi(r,t) is separable as (j)(r,t) — 
g(r)f(t). This leads to the following equations for g(r) 
and f(t), respectively: 



ihv F tr ■ Vg(r) + V l g(r) = Eg(r), 



(4) 



where k l m = sgn(E l m ) (E^) 2 — (k y ) 2 is the wavevector 
of the n-th sideband in the a;-direction, ip l is the phase 
of the oscillating potential in the i-th region. Note that, 
here E l m = (E F — V l + mtuj) / (hv p) has the same dimen- 
sion as the wavevectors k % m and k y . According to Eq. [71 
the Floquet energy Ep can be determined up to an arbi- 
trary integer multiplied by fwj. We can define Eq = Ef, 
E n = Ep + nhtu. By combining the the solutions for f(t) 
and g(r), we can write the solution for the Floquet state 
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Note that, in the free region, i.e., V{ = 0, J n _ m (0) = 
S myn and Eq. ((HJ can also describe the electronic states 
correctly. Here we have obtained the wave functions that 
valid both in the free and oscillating regions in a unified 
form. Then we can use Floquet scattering-matrix formal- 
ism (see Appendix. [X]) to calculate the pumped current. 



B. Pumped current 

Using the distribution function fa Ut (E) for outgoing 
particles and f™(E) for incoming ones, we can find the 
directed current I a in the lead a: 
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FIG. 1: Schematic diagram of a series of time-periodic potentials. 



where, the can be calculated by using the S- 

matrix Sf that we have already obtained in Appendix. El 

fT\E) = Y. E \SF, a p{E,E n )\*ff{E n ). (10) 

f3 pro—n 

The subscript pro-n indicates that the second summation 
in Eq. (ITU1) only includes the propagating states with real 
wave vectors that contribute to the current. The bound 
states near the oscillating scatterer are neglected for elec- 
trons in these states are confined around the boundaries 
without transmitting far away. Substituting Eq. (|10[) into 
Eq. ([9]), we have: 

t poo fE/hvp 

Ia= 6 -^ / dEdk y Y: E 

ntl J-°oJ-E/hv F pro _ n 

\S F ME,En)\ 2 (f}T(En) r:\E)). (11) 

The pumped current at the low temperature limit is: 

Ia = ^ / dEdkyY, E 

nn JE Fermi -Nnw J-E/hvp p pro-n 

\S F ME,E n )\ 2 (flT(E n ) - f™(E)), (12) 

where, a is the outgoing terminal, /3 is the incoming ter- 
minal, N is the maximum index of the sidebands that 
need to be included in the calculation and it is deter- 
mined by the amplitude and frequency of the oscillation 
(proportional to Vi/hia). Next we consider the adia- 
batic limitation: uj — ¥ 0. In this situation, the Floquet 
scattering matrix is almost energy independent when E 
changes by NTiuj. We can rewrite the Floquet scattering 
matrix in adiabatic approximation as SF,a/3(E, E n ) — > 
SQ, a p,n{E), and the difference of the Fermi distribution as 

f(E n ) - f(E) ^nhuj f[E ^Z^ E) -¥ nUuj^§^-, where n 
denotes the n-th sideband in which the electrons are scat- 
tered to. Inserting the above expressions into Eq. (ITT1) . 
the pumped current in the adiabatic pumping regime can 



be reduced to, 

eL v io f°° rE/hv F 

Z7F J-coJ-E/hvF p pro _ n 

n\S ^ n (E)\ 2 ^^- (13) 

We note en passant that, at the low temperature limit, 
only the particles close to the Fermi level can contribute 
to the adiabatic pumped current as indicated by the 
derivative of the Fermi distribution in Eq. (|13j) . We can 
also use the scattering matrix to calculate some other 
important physical quantities such as the directed heat 
current, shot noise, Wigner delay time. 0, MEM 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

A. Transmission asymmetry 

In this section, we use the above formalism to investi- 
gate quantum tunneling through double oscillating bar- 
riers schematically illustrated in Fig. [T] with only the first 
two oscillating barriers included. It is clear that a single 
time-periodic barrier with time-reversal symmetry can 
not produce a pumped current. Here we consider only 
two oscillating potentials which are sufficient to present 
a clear physical picture. The heights of the two barriers 
oscillate in time with the same frequency uj, amplitude 
V\, but with a phase lag A.ip. Thus for a pump with 
Aip ^ nn (n is a integer), the time-reversal symmetry is 
broken and a nonzero pumped current can be produced. 
To understand the origin of the directed current when 
the chemical potentials at both sides of the scatterers 
are equal, we first examine the transmission asymmetry 
for incoming waves in a certain mode Eq — E from both 
sides of the scattering region. 
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T net = J2 {\S F MEn,E)\ 2 ^rf(E) 
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-\S F , 0a (En,E)\ 2 ^fT(E)}, (14) 







where a denotes the right terminal and (3 denotes the left 
terminal. 
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FIG. 2: The transmission asymmetry as a function of the 
phase lag Aip at different incident angles and energies, with 
Vo = 200 meV, Vi = 9 meV, hco = 3 meV, h = h = 50 nm 
h = 25 nm 

In Fig. [2j we plot the transmission asymmetry prob- 
ability as a function of the phase lag Atp at different 
incident angles and energies. Anisotropic tunneling be- 
havior is observed in this time-periodic case. Note that 
for normal incidence, T net is always zero. This phenom- 
ena is the result of Klein tunneling i.e., the incident elec- 
trons at both sides of the scatterer can perfectly trans- 
mit to the opposite side resulting in the vanish of the 
transmission asymmetry. In Fig. [3Ja)-(c), we present the 
transmission asymmetry probability as a function of in- 
cident energy under a certain Fermi energy Ep = 0.2 
eV, while the widths of the barriers are different. Our 
calculation shows that T net vanishes gradually as the in- 
cident energy decreases, since the influence of the oscil- 
lating potentials that break the time-reversal symmetry 
become negligible when the incident energy is far bellow 
Vo as well as the oscillating potentials. The transmission 
asymmetry probability oscillates more seriously when the 
barrier is thicker by comparing Fig. GHa)-(c). This is re- 
lated to the interference between the multireflections of 
the electron waves inside the barriers. So it shows more 
pronounced resonant tunneling behaviors as the barrier 
width increases. We also investigate the effect of the os- 
cillating amplitude V\ and frequency ui on T net as shown 
in Fig. [U[d)-(f). We find that the transmission asymme- 
try probability can be enlarged both by increasing V\ or 
lj for small values of ui. Such results are obvious since 
larger Vi gives rise to a greater number of sidebands serv- 
ing as the transmission modes, while larger ui gives rise 
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FIG. 3: (a)-(c) The transmission asymmetry T nct as a func- 
tion of incident energy with different barrier widths. Vo = 200 
meV, Vi = 9 meV, huj — 3 meV. (d)-(f) The transmission 
asymmetry T ne t as a function of incident energy with dif- 
ferent amplitudes and frequency of the oscillating potentials. 
h = I3 = 50 nm, Vo = 200 meV. In all panels the incident 
angle and the distance between the two barriers are fixed, i.e., 
</> = 30°, h = 25 nm. 



to larger modulation quantum of energy that an electron 
exchanged with the oscillating potentials during the tun- 
neling processes. In Fig.[3J we find that the transmission 
asymmetry vanishes when the incident energy approaches 
V - This is because the ^-component of the wave vector 
k m of the central band and first few sidebands in the 
barrier are imaginary in the incident energy range closed 
to Vo- Imaginary wave vectors corresponds to evanescent 
modes, that decay exponentially and thus lead to total 
reflection at both sides of the scatterer. [32| The breaking 
of the time-reversal symmetry by the interference of dif- 
ferent traversal modes (i.e., central band and sidebands) 
is attributed to the phase lag between the two oscillating 
potentials. 



B. Adiabatic Pumped Current 

Nonzero transmission asymmetry can give rise to a fi- 
nite charge current. Next we proceed to calculate the 
pumped current in adiabatic limit: us — > 0. The expres- 
sion of adiabatic pumped current in Eq. (I13[) is valid when 
the Floquet scattering matrix is energy independent on 
the scale of the order of Nfrui. The lager the number N 
of the sidebands included, the smaller the frequencies for 
which the adiabatic approximation is valid. As discussed 
above, N is determined by the ratio of the oscillation 
amplitude to the frequency. In this subsection, we con- 
sider a weak pump with N = 3. For a strong pump, N 
is much larger and the adiabatic approximation is valid 
only at smaller frequencies. In Fig. El we can see that the 
transmission asymmetry shows a resonance-like behavior 
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FIG. 4: (a) Energy dependence of the adiabatic pumped cur- 
rent, with A<p = tt/2, to = 0.03 meV, Vi = 3w, Vo = 0.2 
eV h — 60 nm, h — h — 100 nm, and Io, a <i = eL y u) / '(2tt 2 ). 
(b)-(c) The adiabatic pumped currents versus the barrier sep- 
aration and barrier width for fixed Fermi energy Ef = 80 
meV. 



of the order of the width SE ~ 0.02 eV resonance, indi- 
cating that the Floquet scattering matrix changes signif- 
icantly over this energy scale. Therefore we can choose 
an oscillating frequency of to = 0.03 meV, for which the 
adiabatic approximation requirement Nhuj <C SE is sat- 
isfied. Note that at zero temperature the pumped cur- 
rent is dominated by the electrons at the Fermi energy. 
The pumped current Ip_ a d shows similar resonance-like 
behaviors as transmission asymmetry as expected (see 
Fig. Ufa)). Ip.ad oscillates as a function of the Fermi en- 
ergy Ef, becomes smooth and eventually vanishes when 
the fermi energy approaches the static potential height 
Vo- It is much different from that of a conventional quan- 
tum pump based on usual nano-scale semiconductors. 
This remarkable feature can be attributed to the same 
reason as the vanished transmission asymmetry when Ep 
is close to Vq as we examined in the last subsection, i.e., 
the transmission probability in graphene through a high 
barrier is lager than that of a low barrier which is closed 
to the incident energy. Another difference between the 
quantum pumps in graphene and that in the conventional 
two-dimensional electron gas system is that the pumped 
current in graphene can change from positive to nega- 
tive or reversely as the energy increases for a constant 
phase lag, while the direction of conventional pumped 
current is determined by the phase lag of the to oscil- 
lating potentials. We also investigated the influence 
of the geometrical sizes of the oscillating double barri- 
ers. In Figs. IDJb) and (c), we plot the dependence of 
adiabatic pumped currents both on the barrier separa- 
tion I2 and width £1(3)- The adiabatic pumped current 
also oscillates when the barrier separation or width are 
changed. This is caused by the Fabry-Perot resonant 



modes formed between the interfaces due to the multi- 
ple reflections. ll[ These results agree with the previous 
work on adiabatic pumping in graphene, e.g. see (l9j . 
With increasing pumped frequency, the deviation of the 
actual pumped current calculated by Eq. (|12[) from the 
adiabatic one calculated by Eq.(|13p becomes significant, 
and thus the adiabatic approximation of Eq. (|13[) as well 
as the parametric scattering approach used in Ref.pjj 
are not valid anymore. 



C. Non-adiabatic Pumped Current 

Next, we analyze the behavior of non-adiabatic 
pumped current Ip as a function of phase lag Aip by 
using Eq. (|12l) . The incident electrons are homogeneous 
in direction under the same Fermi energy Ep from both 
sides of the scatterer. Fig. [SJa) shows the general char- 
acteristics of quantum pumping in graphene. The oscil- 
lating frequency is set to u> — 3 meV, and thus Nhco is 
of the same order of SE. The adiabatic approximation 
discussed in last subsection is no longer valid. The sizes 
of the double oscillating barriers are fixed for I2 = 25 nm, 
h — I3 — 50 nm in the calculation of this subsection. The 
numerical results show that the pumped current Ip is si- 
nusoidal or cosinoidal dependence on phase lag Aip. We 
find that the static barrier V l can also effectively control 
the amplitude and direction of the current as shown in 
Fig. [5jb). The height of the static barrier V l plays an 
important role in determining the pumping, since it is 
related to the component of the wave vector along the 
x-direction in the barrier region. For electron in the n-th 
energy sideband \E+nhuj — V l \ < hvpk y , the correspond- 
ing x-component of the wave vector k % n becomes imagi- 
nary. Such electrons are in evanescent sidebands and 
do not directly contribute to the current. If the static 
barrier becomes higher, the electronic states outside the 
barrier match the hole states inside the barrier. Then the 
transmission modes in the barrier come back to propa- 
gating states, resulting in an increasing overall pumped 
current amplitude. This phenomena in graphene is dif- 
ferent from the transmission in ordinary semiconductor 
two dimensional electron gas, while there is no available 
hole state inside the barrier contributing to ordinary tun- 
neling. One can see in Fig.[SJb) that the pumped current 
oscillates when the static barrier potential is increased. 
This is caused by the Fabry-Perot resonant modes formed 
in the double barrier structure due to the multiple reflec- 
tions, [llj 

The pumped current can also be tuned by the oscil- 
lating amplitude Vi and frequency oj as shown in Fig. [6] 
We find that the pumped current increases with the os- 
cillating amplitude V\. This is because a larger V\ can 
transfer electrons to higher sidebands and thus increase 
the transmission modes of which the interferences become 
more intensive. This behavior is much different from that 
of an adiabatic pump. In the limit of w — > 0, the adi- 
abatic pumped current is irrespective of the amplitude 
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FIG. 5: (a) The pumped currents versus the phase lag Aip 
for different Vo, with Ef = 80 mev, uj = 3 meV, Vi = 9 meV, 
Jo = eL y /(10 5 Tvh). (b) The pumped currents versus Vo, with 
Aip = ir/2. 
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FIG. 6: (a) The pumped currents versus the phase lag Aip for 
different Vi, with E F = 80 mev, uj = 3 meV, V = 200 meV, 
To = eL y /(W 5 Tvh). (b) The pumped currents versus uj, with 
Aip = tt/2. 



of the oscillating potentials V\ , however the criterion of 
adiabatic approximation depends strongly on V\ . [lijj 

In Fig. |6jb), we inspect the influence of the frequency 
uj for several fixed amplitude V\. The pumped current 
shows a trend of rise first, then fall after uj exceeding 
V\. For small uj, electrons can gain a larger modulation 
energy quantum from the oscillating potentials with in- 
creasing uj. The interference of sidebands with higher en- 
ergies gives rise to a lager pumped current. With increas- 
ing pumped frequency uj under a constant oscillating po- 
tential Vi, the oscillating double barriers will transform 
from a strong pump to a weak pump. In the case that uj 
is lager than V\, the oscillating pump is too weak to lift 
the electrons from the central band to any sidebands and 
thus leads to the drop of the pumped current. This fea- 
ture dose not exist in adiabatic pumping regimes, since 
the pumped frequency w is always much smaller than 
the oscillating potential V± , and the adiabatic current al- 
ways increases monotonely with uj. [I?) Therefore, the 
non-adiabatic pumped current can be tuned by changing 
the phase lag, barrier height, oscillating amplitude and 
frequency. 



IV. CONCLUSION 

In summary, we applied the iterative Floquet 
scattering-matrix method to the Dirac system, which is 
able to calculate the pumped current induced by an arbi- 
trary number of parametric oscillating potentials between 
two electron reservoirs without external bias. Based on 
this approach, we investigate theoretically the quantum 
pumping driven by a series of time-periodic potentials in 
graphene monolayer both in adiabatic and non-adiabatic 
regimes. We have determined how the pumped current 
depends on the phase lag, height of the barrier, frequency 
and amplitude of the oscillation. Finally, the physical 
mechanism of quantum pumping is attributed to the 
interference of different central band and sidebands of 
which the phases factor break the time reversal symme- 
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Appendix A: Transfer-matrix and Scattering-matrix 

In Sec. Ill A[ we have obtained the wave functions for 
all the free and oscillating regions. The wave functions 
(see Eq. (jH))) must be continuous at each boundary x — 
x a ' ' "'. At the boundary between region i and i + 1 
(x = 2g +1 ), the coefficients matrices can be related by 
the transfer matrix M(i, i + 1), 

g;)=M(M + l)(^), (Al) 

where B (C) is the probability amplitudes column vector 
of the right (left)-going waves (see Fig. [T]). Matching the 
wavefunctions at £q +1 and realizing that e - l ( E o+nHu>)t/h 
with different n are orthogonal, we have: 

n m 

^ EE J »-C)' e "' ( "" ray+I S 1+c ™ 1 ]' (A 2 ) 

n 77i 
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i 



E 1 



ri+l 



n m 
ui+1 



1 ~\ _ ^ — i(n — rn)ip 1 '^ 1 

ft™ t-n,' 



i+l 



(A3) 



where, Z* 



J+ 1 



-cCq, is the distance of the i-th region. In 



order to convert the above two equations into the matrix 
form, we define: 



= e m -Or, 



(A4) 



(P )nm — Jn— m(~Z ) ' e ^ ^ ! (^5) 
(Qk)nm ^ Jn-m(^) ■ e^"""^ • (A6) 



(Ol)r 



,J n—m \ a. J 



Using the boundary conditions Eq. 
we can obtain the transfer matrix: 



_K L ik 1) (A7) 

A2|) and Eq. (jA3]) 



M(*,* + l) 



7 l 
(7 1 )" 1 

pi+l pi+l 

% +1 % +1 



^pi pi 

Qr Ql 



(A8) 



If we define: 

T(i,i+1) = 



pi+l pi+1 



-li+l 



(A9) 



(A10) 



pi pi 

Qr, Q\ J V Q l R L Ql 

we have all the four sub-transfer matrix: 

M n (i,i+1) = (T^^Tn^i+l), 
M 12 (i,i+1) = (7 < )- 1 T 12 (i,i + l), 
M 2 i(*,* + 1) = ( 7 j )T 21 (m + 1), 
M 22 (m + 1) - (7*)r 22 (M + 1). 

So far, we have obtained the transfer matrices of each 
region. But for an evanescent mode, where k l m becomes 
imaginary, the ( 7 *) mm and(7 4 )~J„ will , respectively, de- 
cay or grow exponentially. An approatch to avoid this 
problem is using the scattering matrix formalism. In 
the following, we obtain the scattering matrix for the 
graphene system following Xu's work on quantum anti- 
dot arrays. [3l| 

The scattering matrix S{L, R), connecting the incom- 
ing and outgoing channels, is defined as follows: 



(B R 



= S{L,R) 



C R 



S n (L,R) S 12 {L,R) 
S 21 (L,R) S 22 {L,R) 



B L 

C R 



(All) 



Then, we use the iterative scheme, if we already have 
S(L,i): 



(A12) 



with the help of the transfer matrix (see Eq. (IAip ). we 
can obtain the scattering matrix of the next boundary 
S(L,i + l): 



S u (L,i+l) = [l-Af 1 - 1 1 (*,i + l)S 1 2(i,t)Af 2 i(i,i+l)]- 1 
■M^(i,i + l)S n (L,i) 

Su(L,i + l) = [1- M n 1 {i 7 i + l)S 12 (L 7 i)M 21 (i 7 i + l)}- 1 
■M n \i, i + 1)S 12 (L, i)M 22 (i, i + 1) 
-M n \i,i + l)M 12 (i,i + l)} 

5 12 (L,i + l) = S 2 2(L,i)M 2 i(i,i+l)5 u (L,i + l) 
+S 21 (L,i) 

S 22 (L,i + l) = S 22 (L,i)M 21 (i,i+l)S 12 (L,i+l) 

+S 22 (L,i)-M 22 (i,i + l). (A13) 

We have obtained the iterative relations of the scatter- 
matrices in Eq. (|A13|) . It is easy to write the initial 
S-matrix: Sl.l = !■ Then we can get every desired S- 
matrix Sl,r using the iterative scheme step by step. Note 
that the current associated with a scattered wave is pro- 
portional to the square of the wavefunction multiplied by 
the velocity. In order to ensure current conservation The 
unitariy Floquet scattering matrix by incorporating the 
velocities of different modes in the incident and outgoing 
terminals, 



v 

SF, a p(in,n) = Sap(m,n)\ -j, (A14) 



where a and j3 denote the terminals, m and n label the 
transmission modes. In this research of a series of time- 
periodic potentials, the velocities in the propagating di- 
rection is calculated as follows: 



hdkS, 



E a 



(A15) 



This iterative scattering-matrix method has the advan- 
tage that we can effectively investigate a complicate sys- 
tem consisting of an arbitrary number of parametric 
pumping regions even though evanescent modes may ex- 
ist in such regions. 
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